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^ Abstract 

By means of certain limit technique, two kinds of generalized Darboux transformations are 
' constructed for the derivative nonlinear Schodinger equation (DNLS). These transformations are 

^ ' shown to lead to two solution formulas for DNLS in terms of determinants. As applications, several 

different types of high order solutions are calculated for this equation. 
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1 Introduction 

The derivative nonlinear Schodinger equation (DNLS) [22l [23] 

lUt + Ur^r, + [{\u\'^U)x = 0, (1) 

has many physical applications, especially in space plasma physics and nonlinear optics. It well 
describes small-amplitude nonlinear Alfven waves in a low-/3 plasma, propagating strictly parallel or 
at a small angle to the ambient magnetic field. It was shown that the DNLS also models large- 
amplitude magnetohydrodynamic (MHD) waves in a high-/3 plasma propagating at an arbitrary angle 
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to the ambient magnetic field. In nonlinear optics, the modified nonlinear Schrodinger equation [5], 
which is gauge equivalent to DNLS, arises in the theory of ultrashort femtosecond nonlinear pulses in 
optical fibres, when the spectral width of the pulses becomes comparable with the carrier frequency 
and the effect of self-steepening of the pulse should be taken into account. 

High order solitons describe the interaction between N solitons of equal amplitude but having a 
particular chirp |10| . In the terminology of inverse scattering transformation (1ST), they correspond 
to multiple-pole solitons. In the case of the Korteweg-de Vries equation (KdV), the poles must be 
simple, that is the reason why high order nonsingular solitons do not exist. Indeed we could obtain 
multiple-pole solutions by Darboux transformation (DT), such as positon solutions [20]. The high 
order solitons for nonlinear Schrodinger equation had been studied by many authors \10 \ \28\ [30]. To 
the best of our knowledge, the high order solitons of DNLS have never been reported. The aim of the 
present paper is to show that such solutions may be obtained by generalized Darboux transformations 
(gDT). 

Recently, the rogue wave phenomenon [15] . which "appears from nowhere and disappears without 
a trace", has been a subject of extensive study. Those waves, also known as freak, monster or giant 
waves, is characterized with large amplitudes and often appear on the sea surface. One of the possible 
ways to explain the rogue waves is the rogue wave solution and modulation instability and there are 
a series of works done by Akhmediev's group [TJ O [3]. Different approaches have been proposed to 
construct the generalized rogue wave solutions of nonlinear Schodinger equation (NLS), for example, 
the algebro-geometric method is adopted by Matveev et al [8l[9], Ohta and Yang work in the framework 
of Hirota bilinear method while the present authors use the gDT as a tool [TI]. In the 1ST terminology, 
the high rogue wave corresponds to multiple-pole solution at the branch spectral parameter in the non- 
vanishing background [16]. The first order rogue wave for DNLS was obtained by Xu and coworkers 
recently [26]. However, the high order rogue wave had never been studied. We will tackle this problem 
by constructing generalized Darboux transformations. 

The Darboux transformation [71 [21], which does not need to do the inverse spectral analysis, 
provides a direct way to solve the Lax pair equations algebraically. However, there is a defect that 
classical DT cannot be iterated at the same spectral parameter. This defect makes it impossible to 
construct the high order rogue wave solutions by DT directly. Thus we must modify the DT method. 
In this work, we extend the DT by the limit technique, so that it can be iterated at the same spectral 
parameter. The modified transformation is referred as generalized Darboux transformation (gDT). 

The inverse scattering method was used to study DNLS with vanishing background (VBC) and 
non- vanishing background (NVBC) [H [5] [6] [T^ I17| . The A^-bright soliton formula for DNLS was 
established by Nakamura and Chen by the Hirota bilinear method [24j and the DT for DNLS was 
constructed by Imai [12] and Steudel |27j (see also j26| ) . One key aim of this work is the construction 
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of the gDT and based on it, the high order sohton solutions and rogue wave are obtained. In addition 
to above two kinds of sohitions, new A^-sohtons and high order rational solutions are also found. 

The organization of this paper is as follows. In section 2, we provide a rigorous proof for elementary 
DT for Kaup-Newell (KN) system. Furthermore, based on the elementary DT, we construct for KN 
system the binary DT, which is referred as the DT-II while the elementary DT is referred as DT-I. We 
also iterate these DT and work out the A^-fold DT's both for DT-I and DT-II, and consider two different 
kinds of reductions of the DT of the KN system to the DNLS. In section 3, the generalized DT-I and 
DT-II are constructed in detail by the limit technique. In section 4, we consider the applications of 
the generalized DT (gDT) and calculate various high order solutions, which include high order bright 
solitons with the VBC, and high order rogue wave solutions. Final section concludes the paper and 
offers some discussions. 



2 Darboux transformation for DNLS 

Let us start with the following system — Kaup-Newell system [13 

m + Uxx - = 0, 



(2) 



-\vt + Vxx + = 0, 

which may be written as the compatibility condition 

Ut-Vx + [U,V] = 0, (3) 

of the linear system or Lax pair [18j 

$x = f/^, (4a) 
$t = (4b) 

where 

U = [-^as + ^Q], V = -^as + - ^<3^^3 + ^Q^ - ^Qxcr^, 

with 

These equations ([2]) are reduced to the DNLS Q ioic v = —u*. For convenience, we introduce the 
following conjugate linear system for Q, 

-^x = ^U, (5a) 
-qi^ = ^^)V. (5b) 
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2.1 DT-I 



In the following, we first consider the DT of the unreduced linear system (jl]). Generally speaking, DT 
is a special gauge transformation which keeps the form of Lax pair equation invariant. The explicit 
steps for constructing DT in 1 + 1 dimensional integrable system are as following: first, we consider 
the the gauge transformation 

where Di and are unknown matrices which do not which depends on C,. Then imposing that ^[l] 
is a DT we have 

D[l], + D[l]U = U[l]D[ll (det(Z)[l])), = Tr(f/[1] - U) det(Z)[l]), 

where U[V\ represents transformed U matrix. After some analysis, we find the following elementary 
DT (eDT) for g]): 

where $i = is a special solution for linear system dH) at = Ci, and ^*^cti is a special 

solution for conjugate linear system ([5]) at ^ = — Ci- Here we point out that this DT for DNLS was 
first derived by Imai |12j . With the help of the DT, Steudel[27] and Xu et al [26j calculated various 
solutions for DNLS. Next, we give a rigorous proof that the above transformation does qualify as a 
DT. 

Theorem 1 With $i and D[l] defined above, ^>[1] solves 

$[1], = c/[i]$[i], mt = v\i]^i] 

where 

U[l] = --^^3 + -0[l], Q[l] = cjiQai - 2CicTiPi,x-^i, 

and 

Namely, the D[l] qualifies as a Darboux matrix. Correspondingly, D[l]~^ is a Darboux matrix for the 
conjugate Lax system. 

Proof: To begin with, we notice 

C + u C - u 
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What we need to do is to verify 



U[l] = D[ll,D[l]-^+D[l]UD[l]-\ 
V[l] = D[lltD[l]-^ + D[l]VD[l]-\ 



(7) 
(8) 



First we consider The residue for function Fi(C) = D[l]^^D[l]-^ + D[l]UD[l]-'^ - U[l] at 
C = Ci is 



Resc,(Fi(C)) = 2Ciai [-(/ - + (/ - Pi)C/(Ci)Pi] ai = -2Ciai(/ - Pi)$i 



cji = 0. 



where we used the relation = — Similarly, the residue of function -Fi(C) 

at C = — Ci is 

Res„Ci(i^i(C)) = -2Cif7i[-Pi,,(/-Pi)+PiC/(-Ci)(/-Pi)]ai = 2CifTi 
Due to 



ai^i{I-Pi)ai = 0. 



1 



1 



and 



[/[!] = -^a3 + -Q[l], 



Q[l] = aiQcTi - 2CicriPi,a;C7-i, 



(9) 



the function Fi{()D[l] is equal to zero at C = 0. Thus the function -Fi(C) is analytic at ^ = 0. It is 
easy to see that -Fi(C) — )■ at ^ — oo. Therefore the equality ([7]) is valid. 



Now we turn to the time evolution part ([5]). We introduce a matrix V[l] 



C73 + ^Q[l] + 



+ ^Vi, so that F2(C) = D[l]^tD[l]~^ + Z?[l]FL>[l]-i - V[l]. Proceeding similarly as above, it is 
found that F2{C) is analytic at C = iCijO ^iid tends to zero as C — )■ oo, thus F2{C) = 0. 

In the following, we show V[l] = V[l]. Because of the compatibility condition (Z)[l]$)a;f = 
{D[l]^)tx, we have 

u[i]t-m,+ [u[i],m] =0. (10) 

Identifying terms of 0(C) in pO|) . we have 

W3,V2]=0, (11) 

[as,Vi] = [V2,Q[l]] + 2Q[lU (12) 
[Qm,Vi] = V2,.. (13) 
From (jlip. we have V2^^ = 0, where V2^^ denotes the the off-diagonal part of V2- Similarly, we have 

V,'>^ = iasQ[l]x-^as[Q[l],V2] (14) 
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through ([11]). Substituting into (fT3|) and solving it yields 

V2 = -ig[i]V3 + f{t). 

Letting Ci = 0, one can readily obtain V2 = —laiQ^aia^. Therefore we have f{t) = 0. Moreover, 
Vi = Q[lf - \Q[l]^(j^. 

Similar argument could show that does qualify as a Darboux matrix for the conjugate Lax 

system. Thus the proof is completed. 



Remark 1 In addition to Q there is a different representation for 

Q[l] = ai{I - 2Pi)Q{I - 2Pi)ai + ^a^a^il - 2Pi)ai. 



(15) 



which is appeared in papers \21\ \26^ . We will use Q rather than ()15p . since the former is more 
compact. 



To derive the A^-fold DT for this elementary DT which is referred as DT-I, we rewrite it as 

D[l] = 

Assuming different solutions $j = ((/>i, ^pi)'^ of ([4]) at (" = Ci a-re given, we may have 



Proposition 1 21.\M^ The N-fold DT for DT-I can he represented as 



D^ = D[N]D[N -!]■■■ D[l] = C^^a^ + h ° ] a^C^ 

fc=0 \ PkJ 

where Ui are determined by the following equations 

aov?^ + aiCi</'i + • • • + a2i_iCf "Vi + a2zC^Vi = "C^'^Vi, when N = 21 + 1; 

aofi + aiCi(t)i H h a2i-2Cf~'^^i + a2;-iC^'~Vj = -(1 vt, 

i = 1,2,- ■■ ,N. And Pi = Oiiipj o (j = 1, 2, • • • , A). 



^2'-^- when N = 21. 



The transformation between the fields is the following: 
1. When N = 21 + 1 

"det(S)" 



u[N] = -V 



det(A) 



, v[N] = -u + 



where A = [Af , A^,-- - ,Al),B= {BJ , BJ, ■■■ ,B%), 

A, = {^,Xi(Pi,--- ,cf"Vi,cfV*), 



(16) 



(17) 



Bi = {^,Xi^i,--- ,cf-V^,cf+Vi)- 



2. When N = 21 



u[N] = u ■ 



detjD) 
det{C) 



v[N] = v + 



det{D{ipj ^ (pj)) 



where C = {Cf, , ■ ■ ■ ,C]^) , D = {dJ , D^,--- ,D 



2 5 J^NJ> 

2/-2_ a21-1. 



(18) 



2.2 DT-II 

In this subsection, we will show that the so-called dressing-Backlund transformation \18\ 129]. denoted 
by DT-II in the present paper, may be constructed from above DT-I. For convenience, we rewrite D[l] 
and L'[l]^^ as following 



D[l] 




1 



c 



(19) 



Suppose another solution = (xijV'i) for the conjugate system ([5]) at C = ?i is given, then ^'i[l] = 
^iZ)[l]~-'^|(^=g^ is a new solution for conjugate system (^'[1], U[l\^ ^[1]) at ^ = ^i. It is easy to see that 
cri^'i[l]-^ is a special solution for Lax pair ($[1], f^[l], ^[1]) at C = — .^i. Therefore, we could construct 
the second step DT D[2\ by the seed solution cti^i[1]"^. By direct calculations, removing the factor 
C -il, we have the DT-II 



T[l] = I + 



A 



c - 6 c + er 



A 



2 lo /3, 



1^1, 



(20) 



where 



a 



^1 



6 

.0 Ci 



.0 ii, 



$1. 



Furthermore, we have 



r[i]-i = / + 



B 



B 



C-Ci C + Ci 

The transformation between Q and Q[2] is 

Q[2] = Q + [A - a^Aasl 




(21) 



(22) 



Above discussion indicates that the DT-II T[l] is indeed a two- fold DT for D[l] in the case of DNLS. 
We remark that in the case of two component DNLS the analogy of DT-II exists [19] while the 
corresponding DT-I has not been constructed. 
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In what follows, we consider the iteration for the DT-II T[l]. Assume we have N distinct solutions 
^ii^J'i) = V'i)^ of dH) at C = and N distinct solutions ^'j(fj) = (Xi^'^i) of ([5]) at C = Mj. Similar 
to DT-I, we work with DT-II T[l] and have the following proposition 



Proposition 2 The N-fold DT for the DT-II could be written as the following form 

N 



Tn = T[N]T[N -!]••• r[l] =1 + ^ 



Ci (TaCjfJa 



and 



T^i = T[l]-^T[2]-^ ■ ■ ■ T[N]~^ = ^ + X] 



N 



(23) 



(24) 



Proof: We calculate the residues for both sides of 

An 



Res 1^=1., (Tat) 
Res|^=_i,^(rAr) 



{1 + 



Ui - UN Ui + VN 



-{1 + 



A 



N 



) • • • a^AiG^ ■■■{I 



0"3^lO-3 

Vi -Vi Vi + Vi 



£73^1 0-3 



)• 



-Ui - UN -Vi + Vn -^i - ^1 -^i + Z^l ' 

Because of Res|^=^. (T/v) = — o"3Res|^=_,y. (Tjv)cj3, equation (p3]) is valid. Similarly, (p^ can be proved. 

The N-fold DT-II Tn allows us to find the transformations between the fields n[0], v[0] and u[N]^ 
v[N]^ which are given below 



Theorem 2 The N-fold DT-II Tn induces the following transformations for the fields 



u[N] = m[0] - 2 



detMi 



detM 

where M = {Mij)NxN, N = {Nij)NxN, Mij -- 
( 



v{N\ = v[0] 2 



detiVi 
detiV 



(25) 



Ml 



Mil Mi2 
M21 M22 



Mi^ Vi 

M2N ■02 



\ / iVii Ari2 



Mni Mn2 ••• Mnn IpN 
\ <Pi 4>2 ■■■ 4>N J 

Proof: Since Tat ([23]) is the N-fold DT of (HD, we have 



iV21 N22 

Nni Nn2 
\ fl ^2 



■ ■ ■ NiN Xl ^ 
• • • N2N X2 

■ ■ ■ Nnn Xn 



Tn,. + TnU = U[N]Tn. 



It follows that 



N 



Q[N]=Q + Y,[C^-cJ3Cia3l. 



i=l 
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Thus, we need to calculate the explicit forms for Cj. Proposition 2 implies Cj = Res|^=^. (T/v), it 
implies that Cj's are the matrices of rank one. Thus we may assume Q = \xi){yi\. Similarly we may 
set A = \wi){vi\. 

On the one hand, because of TjyT^^ = I, we have 



(26) 



where the fact that the residue T^Tj^^ at C = equals to zero is taken account of. On the other 
hand, we have 

Noticing that the rank of T^^\(^=^- equals to 1, we may obtain 



{yi\ = *z 



Now substituting {yi\ into ([26]) leads to 

N 



i=l 



fJ-i + i^i 



0, {l = l,2,...,N). 



(27) 



Solving ([27|) gives us 



[kl),|x2),--- ,|xAr)]i = (/)2,--- at] M \ 

[\xi),\x2),--- ,\xn)]2 = ['Pl,^2,--- ,^n]N-^, 

where subscript i and 2 stand the first and second rows respectively. Finally, the relations between 
the fields can be represented as 



u[N] = u[0] + 2 



v[N] = v[0] - 2 



■02 
X2 

\xnJ 



u[0] - 2 



det Ml 
detM 



v[0] + 2 



detiVi 
detiV 



This completes the proof. 



9 



2.3 Reduction 

So far we have been working with the DTs for the general Lax problem ^ and certain solution 
formulae have been given for the system ([2]). However our main task is to construct solutions for 
DNLS ([T]), therefore we have to consider reduction problem. It is easy to see that two reductions 
V = u* and v = —u* are simply related [27], so we may consider either of them. For the DT-I, let us 
assume v = —u* or = —Q, where ^ denotes the complex conjugation and transpose. To implement 
the reduction, we need to choose the seed solutions properly. Indeed, assuming 

CiGiH, and (^i = (28) 

then Q[l\, defined by Q, satisfies the reduction relation = — The DT ^ with the reduction 
condition (|28p may be employed to construct bright or dark solitons of DNLS with the non-vanishing 
background. 

Let us now turn to the reduction of the DT-II. Assuming Q = and 

ei = Ci*, and {xi,i^i) = i^l,vl), (29) 

([22]) yields Q[2] = Q[2]l 

To iterate the DT-I and DT-II, we must verify that they keep the reduction conditions (j28p and 
(j29p . The latter merely depends on the symmetry of equations (jH, thus it holds automatically. For 
the former ()28p we claim that 

Proposition 3 Both DT D[l] andT[l] keep the reduction condition (|28p invariant. 

Proof: Direct calculations. 

Due to above analysis, both DT-I and DT-II may be reduced to find solutions for DNLS. However, 
the DT-I under ()28p is conveniently used to construct the N-dark or bright soliton solutions of DNLS 
with NVBC, while DT-II with ()29p may be properly adopted to represent the N-bright solitons and 
N-breathers of DNLS. 

3 Generalized Darboux transformations 

In this section, we construct the corresponding generalized Darboux transformations (gDT) associated 
with D[l] and ^[l]. We will follow the approach proposed for the nonlinear Schordinger equation in 
[llj . Indeed, while both DT-I and DT-II considered above are degenerate at = Ci in the sense that 
^[l]lc=Ci*^i = ^[l]lc=Ci*^i = 0, we may work with 
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or 



4>i = ,i„™fi)ba±i 

which serves the seed sohition for doing next step transformation. 



3.1 gDT-I 

To construct the gDT associated with DT-I, we assume that n solutions {ipi,cl)i)'^ are given for the 
Lax pair at C = Ci = li " " " j First, we have the elementary DT 



D 



-Ci— c 



As observed above, by virtue of the limit process, we find that 

M°\=0+^^i /^i(Ci + e)\ ^[Oli d (MO 



lim ■ 



Vl(Cl 



is a non-trivial solution for Lax pair (jl]) with n = n[l] and = at C = Ci) which may lead to the 
next step transformation 

This process may be continued and we have the following theorem. 



Theorem 3 Let {ipi,cf)i)^ be the solutions of Lax pair at C, = C,i {i = 1, 
possesses nii order zeros at C, = Then we have the following gDT-L 



, n) and assume DT-I 
(30) 



where 



and 



i=l 




(31) 



and 



[j-2] 



up • • • M| 



[mi-l] 



[k] 



0-1, 



D 



if 5^ = I; 
if 5f = 0. 
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Proof: To construct the gDT-I, we start with the eDT 



By means of the nontrivial solutions (/>i[l]), we may do the next step of transformation D^^\ 

Taking account of the given seeds {ipi, 4)1)^ , we perform the following limit 



9?^ 



i 



lim 



which yields the formulae presented in above theorem. This completes the proof. 

To have a compact determinantal representation for the gDT-I, we may take the limit directly on 
the N-fold DT-I (fT6]) . It follows from p!7|) and (fT8|) that the transformations between the fields are: 



1. When N = 21 + 1 

u[N] = -v[0] 



det(B) 



det(^) 



v[N] = -ti[0] + 



det ^ 4>j)) 
det(yl((/3j o (j)j)) 



(32) 



where 



A I aT ^ aT dJ^^ 



-A 



T 



... aT ^aT ... " 



-At 



B=[Bl,-B 



d 



jmi — 1 



5 



(mi - l)!dC'"i-^ ^ 



_ . . . 



d 



•vrin — l 



(mn - l)!dC 



and Aj, Bi are the same as (H 
2. When TV = 21 



n[iV] = u[0] 



where 



det(L>) 
det(C) 



det(C(v3j •(-)■ (/>j)) 



(33) 



^"'^^^'^^^^'■■■'(m„-l)!dC 



d 



■mn-1 



m„-l 



(m„ - l)!dC 



(m„-l)!dC'""-i " ' ' 



T 



(m„ - l)!dC""-i " 



and Cj. Di are the same as (fTSl). 



Thus, we complete the construction of gDT-I for ([5]), which could be considered as a generalization 
for DT studied in [l2l [261 [27| . 
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3.2 gDT-II 

In this subsection, we consider the generahzation for DT-II. To this end, we assume that n solutions 
^i(C = /^) = {4'i-,^i)'^ are given for the Lax pair at fi = fii and n solutions ^i{C = v) = (xi^i^i) are 
given for the conjugate Lax pair at v = Vi (i = 1, ■ ■ ■ , n) . 



Theorem 4 Let <I>j be the solutions of Lax pair at C = fJ-i o.nd be the solutions of conjugate Lax 
pair at C, = Ui (i = 1, • • • , n), 

r 
i=l 

assume DT-LI possesses mi order zeros at Q = it/ij and inverse of DT-U possesses mi order zeros at 
Q = ibfj. Then we have the following gDT-II 

where 



>[j] 



C - /ij C + Mi ' 



/3 



a- 



and 



\0 fiij \0 UiJ 



^(j-/)!d^^-'^^"^=^- 



[m.i-l] 



d^- 



A/ 



1=0 



1 

~2 J • 



if 51 = 2; 



Ai= J2 Nt'^...Nl'^...Nl'^^-'^...Nf>\ 

if6{ = 2; 



i^6i = 0. 



7^[i]-i| 



iJ6{=0. 
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Proof: Noting that the DT-II is given (j20p and using the hmit technique, we could obtain special 
solutions for new Lax pair (g]) ($[1], f/[l], ^[1]) at C = Ml and conjugate Lax pair ([5]) (^'[1], y [1]) 
at C = z^i, i.e. 



[1] 



[1] 



Tl^\^=^^+S^i{fii+6) _ [0] d ,q^(„^ 
lim ^ -I^ \(=f,l-^9l\^,=^ll+ bi9i(Pi), 



lim 

(5-5.0 



^i(z.i + (^)rfi \=,,+5 _ d 



where 



oi — 5 9 , ill — — 



1/2 - /if 



Therefore, we could continue to construct the DT-II T[2] for the new system 

^ C - C + z^i ' ^ C - Ml C + Ml ' 

Generally, taking account of the given seeds <I>j and ^j, we perform the following limit 



lim 



lim + 5) 

(5— s-O 



[mi-1] 
1 



• • • T ' 



1 j IC=A'i+'5 



<5J 

[0]-l [mi-l]-l 



. . . t' 



1 . . . 



(5^ 



and mathematical induction leads to gDT-II (j34p . This completes the proof. 
Due to above proposition, the transformations between the fields are 



n rrii — l 



u[N] 



v\m 



i=l 3=0 
n rrii — l 



(35) 
(36) 



i=l j=0 

As before, the above formulas (j35|) and (j36p could be rewritten with the determinant form, i.e. 



u[N] = u[0] - 2 , ^[iV] = ^;[0] - 2 



(37) 



where 



I' p[ll] p[l2] . . . p[lr] 
p[21] p[22] ... p[2r] ^ 



p[rl] p[r2] 

\ 4>i 4>2 



^ 



p 



( p[ll] p[12] 
p[21] p[22\ 



p[lr]\ 
p[2r] 



^p[rl\ p[T2] ... p[rr] ^ 
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with 



and 



with 



5 ; 



{mj-iy.a^rnj-l9j 



P, 



kl 



{k - !)!(/ - 1)! dv^-^dii^~^ V + 

g[21] g[22] ... g[2r] 3^ 



fl — V 



Qi 



gH] Q['^2] ... g[rr] ^ 

y ^ ... ^ 



g[21] g[22] 



Xi — [Xi, d^Xi} ■ ■ ■ ) (rrii-iy. du"H-'^ Xi 



T 



1 a^j- 



Qk+l-2 /^.(^,)^3$.(^) 



fl=fli,U = l'j 



\^g[rl] Q[r2] ... gH^ 



(fc- 1)!(^ - 1)! dv'^-^d^i^-^ 



+ 



fJ.=fli,U=Uj 



+ V — V 

According to above theorems, it is not difficult to see that the reductions (j28p and (j29p are still 
valid for gDT-I and gDT-II. For the gDT-II, to reduce system (jl]) to DNLS, we must set /ij = v*. 



4 High order solutions for DNLS 

Integrable nonlinear partial differential equations are well known for their richness of solutions. To 
construct those solutions, a number of approaches have been proposed including Inverse Scattering 
Transform (1ST), Dressing method, Hirota's bilinear theory and Darboux (Backlund) method, etc. 

While classical DT is known to be a convenient tool to construct N-soliton solutions, it can not be 
directly used to obtain the high order solutions, which correspond to multiple poles of the reflection 
coefficient in the 1ST terminology (see [28] and the references there). We will show in this section that 
the gDT derived above can be applied to obtain various solutions for DNLS. Indeed, apart from the 
high order solutions, a kind of new N-soliton solutions will also appear. In subsection 4.1, we consider 
the solutions with the VBC. N-rational solitons, high order rational solitons and high order solitons 
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are worked out. In subsection 4.2, we construct the high order solutions with NVBC which include 
the high order rational solutions with NVBC and high order rogue wave solutions. 



4.1 Solutions with VBC 

Applying DT-I to vacuum, we may obtain three kinds of solutions, namely plane wave solutions, N- 
phase solutions (periodic solutions) and N-soliton solutions (see |26j). Additionally, if we take limit 
of the soliton solutions, we can find rational solutions |13^ 126] . In this subsection, we consider the 
N-rational solutions first. As we pointed out, the rational solitons are the limit cases to the soliton 
solutions. The different behaviors of the high order rational solitons and high order solitons are 
indicated. 

In the first two cases, gDT-I will be used, while for the case 3, it is more convenient to use the 
gDT-II since the spectral parameters need be conjugated each other. 
Case 1: N-rational solutions. 

We first consider the rational solitons and their higher order analogies. For the seed solution n = 0, 
the special solution for Lax pair Q with the reduction u = —v* is 



where c is a complex constant. For higher order solutions, c will be taken as a polynomial function of 
so that the high order solutions with free parameters may be constructed. 
To obtain the N-rational solitons, we introduce vectors 




(38) 



■2N-2 




Z = 



2N-2 



2N-1 



and define the matrices 





Y = 



Z = 



where yi = y\t^=Q^c=c, and Zi = z\q=q^c=c,, th 
the N-rational soliton can be represented as 





(39) 



Taking Ci = io, we have 



2i(a^3:-2t + a^e) 



u[l] 
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I One soltton-l Two-soliton (1,2) — - One soliton-2| 



Figure 1: Two-soliton (1,2) with the parameters Ci = 2i, ^2 = 4i and ci = C2 
parameters Ci = 4i and ci = 1; One soHton-2 with parameters Ci = 2i and ci 



1; One soliton-l with 
1. 



which appeared already in [26]. The velocity for this rational soliton is a^/4 and the center is along the 
line a^x — 4t -|- c?c = 0. The altitude for |u[l]p is 16/a^. A simple analysis shows that this 2-rational 
soliton does not possess phase shift when t — t- ±00, which is different from the 2-soliton of NLS. This 
phenomenon is illustrated by Fig. [TJ 
Case 2: High order rational solitons 

Next we consider the high order rational solutions. Set the matrices 



y? 



\ 



{2N~1) 
1 



( 



^(2Ar-l) 



where the superscript represents the i-th derivative with respect to It follows that the A^— order 
rational soliton for DNLS with VBC can be formulated as 



detFi 
det Z\ 



(40) 



By choosing appropriate parameters, we have the first and second order rational solitons with VBC, 
which are plotted in Fig. [2j 
Case 3: High order solitons 

To obtain the high order soliton solutions, we start with the seed solution ii[0] = 0. The special 
solution for Lax pair Q with the reduction u = —v* at n = u[0] is 

/g-iM-2(a::+2M-2f+c)\ 
\ iAt-2(x+2/i-2i+d) j ' 
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10 -10 



(a) 1-st rational soliton (b) 2-rd rational soliton 

Figure 2: (Color online): High order rational solution with VBC: The parameters Ci = 2i and ci = 

and the special solution for conjugate Lax pair ^ with the reduction u = —v* at u = u[0] reads as 



{x, 



Then, the N-th order soliton solution for DNLS is given by 

'detMi 



u[N] 



detM 



M = {Mij)NxN 



(41) 



where 



M, 



1/2 — /i^ 



Ml 



d/i 



N- 



-if], Y = [i;,—ip, 



d 



N-l 



dz/"^' 'dl/^ 





The second order and third order soliton solutions are shown in Fig. [3l The high order soliton 
with more free parameters may be obained if c or d are taken as polynomial functions of fi. 



4.2 Solution with NVBC 

The solutions with NVBC may be obtained by applying DT to zero solution. As illustrated in |27j . 
one fold DT-I could be used to yield the plane wave solution. Thus we will consider the high order 
rational soliton solutions resulted from vacuum first. To find more general solutions with NVBC, we 
may apply DT to the general plane wave solution. This will be considered in case 2 and the genuine 
rational solutions and their high order analogies with NVBC will be calculated. In case 3, we construct 
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2 4 




-^10 



-60 -10 




(a) 2-nd order soliton (b) 3-rd order soliton 

Figure 3: (Color online): High order soliton with VBC: The parameters n = 1 + i and c = d = 

high order rogue wave solutions. As above, gDT-I will be employed in first two cases and gDT-II will 
be adopted in the last case. 

To the high order rational solution with NVBC and high order rogue waves, because the all locate 
in the branch point spectral. We must use some tricks to deal with this problem. Besides the high 
order rogue wave solution, the high order breather solution and periodic solution can be readily to 
obtain similarly. Because they are nothing but using the the formula p7p directly like the above 
subsection: case 3. Thus we omit them in our work. 

Case 1 : High order rational solitons with NVBC from vacuum 

To obtain high order rational solutions with NVBC, the order of determinants must be odd. Define 
the matrices 



where 



( ] 























-27V-1. 



and (j), tp are given by ()38p . Then the high order rational solitons with NVBC can be represented as 

/ detyi\ 



u[N] = 

\ det Zi I 

The first and second order rational solitions in NVBC are shown in Fig. [H 



(42) 



19 




I«l 




(a) 1-st order rational soliton 



(b) 2-nd order rational soliton 



Figure 4: (Color online): Rational solitons with NVBC: The parameters (^i = 2i and ci = 

Case 2 : High order rational solitons with NVBC from plane wave solution 
To construct these solutions, we take the seed solution as 



u[0] = Aexp(2i02), 

where 62 = h[ax — {A^a + o?)t -|- c], c G M. The corresponding fundamental solution for Lax pair (jl]) is 



where 



and 



' gi(e'i-H6»2)+0 g-i(6>i-6»2)-<^^ 



1 2 + aC\ 



(43) 



-^/-{2C~'^ + aY -AA^C~'^[x -{a + A^- 2C^)t + d]. 



To resolve the reduction (j28p . we assume 

_4^2^2 _ (2 aC^f > 0, d G M, C G The corresponding 
dark solitons and bright solitons were studied and analyzed in [27\ [26] in details. The authors of [27] 
also illustrated certain limit cases, but they did not give the explicit expression for those solutions. 

In the case —AA'^(^'^ — (2 -|- a^^^s _ q <|) given by the formula (j43p . which does not qualify as the 
fundamental solution, is in fact a constant. Thus the gDT could not generate interesting solutions. 
To obtain meaningful solutions, we must find another solution for ([H). To this end, we turn to the 
limit technique. 
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For convenience, we consider the special case a = c = 0,A = l which leads to the genuine rational 
solutions. We will expand the solution for Lax pair at = i- With the special solution 



$1 



$1 



C=i(l+/) 



c 



fl/2 



c 



at / = 0, we have 



$1 



Yo + Yif + --- + YnP + 
(xn\ 1 d'^ 



lim 



r$i(/). 



/=o n! d/" 

The high order genuine rational solution in NVBC is represented as the following: 
1. When = 2Z — 1, we have 



u[2l - 1] = -1 



det(S) 
det(A) 



(44) 



where 



min(i— l,2jr— 2) j_l 
Ai,2j-1 = Bi^2j~l = i^-'"^ ^ C2j-iXk, (j = 1, 2, • • • ,1 - 1), -Bj,2/-1 = i^' ^ C2iXk, 



k=0 



k=0 



min(i— l,2j — 1) 



^ C2jyk, 



k=0 



n!(m — n)l 



2. When N = 2l,we have 



u[2/] = 1 



det(L>) 
det(C) 



(45) 



where 

Ci,2j-1 



min(i-l,2j-2) 
2j-l \ " rik 



min(i— l,2jr — 1) 



k=0 



Cl.^iXk, a,2j = i'' Yl ^^jVk, (j = 1, 2, • • • , I) 



k=0 



i-l 



L2k 



Q,2fc, (A; = 1,2,--- ,Z-1), A,2/ = i''+'^C2m?/fc- 



fc=0 



In particular, taking d = ef, where e is real number, we have 

/2(2x-6i + i), 



xo 

Xl 

yi 



V2{2x-6t-i), yo 
"2 
3 
"2 



^/2 
V2 



-x^ - Qx'^t + ISxt'^ - 18t^ - 4x + 20t + 2e + i(- - + 6xt - 9t^) 



- 6x^t + ISxt^ - 18t^ - 4x + 20t + 2e + i(x2 6xt + 9*^) 

3 2 
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Thus by means of above formula (|44p . the first order genuine rational soliton solution reads 



n[l] 



{-2x + 6t- i)(-2x + 6t + 3i) 
{-2x + 6t + i)2 



Similarly, (j45p provides us the following second-order genuine rational soliton with NVBC 



u[2] 



L\L2 




where the superscript * stands for complex conjugation and 

Li = -8x^ + 72x^t - 2lQxf + 216t^ + 121^2 - 72\xt + lOSii^ - 18x + 102i + 3i + 12e, 
L2 = -8x^ + 72x^t - 216xt^ + 216t^ + 36ix^ - 216i2;i + 3241*^ + 30x - 42t + 12e - 15i. 

Case 3 : High-order rogue wave solutions 

The rogue wave solution for DNLS was first derived in [26] via DT, to the best of our knowledge. 
However, the classical DT can be not be used directly to obtain high order rogue wave solutions. 
According to above, we can see that the gDT is a very efficient way to obtain high order solutions. 

In order to get the higher-order rogue wave solutions, for simplicity, we consider the seed solution 
u[0] = exp(— ix). The corresponding fundamental-matrix solution for Lax pair is 



and -F(C) is a polynomial function for C,. Like above case, by means of the limit technique, we expand 
the special solution at = 1 -|- i 




where 



a = 



m)-\\ - 2i + K^)f'\ a = exp -\C\x + 2Ch + F{0) , A = (-4 - C^)V2 




at / = 0, and find 



$1 



where 
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Explicitly, we have 
xo = 
yo = 

Xl = 



exp [—-ix]{2x — 2it — 1 — 1), 
exp [^i2;](2x - 2it + 1 + i), 



exp [-^ix\ 



-^x^ + ix'^t + xt^ -^it^ + ^^x^ + {l-i)xt 2 



1. 5 13. 1 „ 1 

--1X X H it t + 2e-\ h 2io 

2 2 2 2 2 ^ 



yi 



exp [hx] 



--x^ + \x^t + xt^ - -it^ 
3 3 

1 



1 + i 



1 + i 
(1 _ i)xt + -—t 



1. 5 13. 1 „ . ^. 

--\x X -\ \t 1 + le h 210 

2 2 2 2 2 ^ 



Therefore, the A^-th order rogue wave can be represented as following 

r ■ 1 /detMi\ 
.[iV]=expH.]-^-^j^, 



where M = (M, 



Ml 



^ Mil Mi2 

M21 M22 

Mm Mn2 

\ Xo Xl 



MiN yo ^ 
M2N yi 



Mnn yw^i 

XN-l 



Mi, 



El ^ \i+j-{k+l+l)r^-k-l 



k=OA=0 



j~{k+l+2) 



\Yla^Yi (1 - 1)^-'=-^ (1 + i)-'-'-^ + iYlYi (1 + ly-'^-^ (1 - iV-'-i 



In particular, the first-order and the second order rogue wave solutions are given by 

[2t^ + 2^2 - 3 - 2i(x + 3t)] [2x2 ^ - 2i(x - t)] 



u[l] 



[2x2 + 2*2 + 1 + 2i(x - 



■ exp[— ix] 



and 



u[2] 



L\L2 



exp[— ix] 
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(a) 




(b) 



Figure 5: (Color online): Second order rogue wave with the parameters (a) e = g = 0; (b) e = and 
g = 100. 



where 



L2 



72 [e^ + 5^] + [48x3 - lUxt"^ + 721^2 + lUlxt - 721*^ - I44x - 72t + 36i]e 

+ [-lUx^t + 48*3 + 72ix2 - 144ixt - 721*^ - 72j; + 432t - 361]^ 

+8x^ + 24x^*2 + 24x2*'' _^ g^6 _^ 24ix^ - 24ix^t + 481x^*2 _ 481x2*3 

+24ixt^ - 24it^ - 12x'' + 48x3* - 216x2*2 + 48xi3 + l80t^ + 481x3 

-2881xt2 - 3361*3 _^ 90^.2 _ j2xt + 666*2 _^ 54-^ _ -^gg-^ _^ 

72[e2 + ^2] + [48x3 - I44x*2 - 721x2 + 4321x* + 721*2 _^ ^^^^ _^ 216* - 1801]e 

+ [-144x2* _^ 4g^3 _^ 2161x2 + 1441x* - 2161*2 + 216x + 144* + 36% 

+8x^ + 24x^*2 + 24x2*4 _^ g^6 _ 24i^5 _ ^2[xH - 481x3*2 - 1441x2*3 - 241x*4 

-721*^ - 60x4 _ 1443.3^ _ 5042:2^2 _ 144^^3 _ gQ^4 _^ 4g-^3 _^ 2881x2* 

+5761x*2 - 5281*3 - 198x2 + 504x* - 486*2 + 901x + 4141* + 45. 



These solutions are plotted In Fig. Owlth different parameters. 



5 Conclusion and Discussion 

The theory of DT Is developed and two generalized Darboux transformations, gDT-I and gDT-II, 
are constructed for DNLS. With the help of them, two generalized determinant solution formulae are 
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obtained for this physically relevant equation. Moreover, high order solitons, high order rogue waves 
and rational solutions are given explicitly. We remark that the gDT-II is still valid for N-component 
DNLS system. In addition, the above formula can be easily modified and applied to so-called Fokas- 
Lenells equation [T7]. 

As shown in Figure 5, the second order rogue waves exhibit dynamics which varies according to 
the different values of the parameters. It is interesting to study the dynamics of the general high order 
solutions. Also, it is interesting to calculate the energy of the high order solutions, which could be 
done directly. However, the calculations involved in is very tedious. These questions may be answered 
efficiently by a detailed analysis of DNLS in the framework of inverse scattering method. 
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